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S_l , We present experimental results on a method to perform polarimetry on ensembles of single photons. Our setup is based on a measurement 

a 

, method known to be optimal for estimating the state of two level systems. The setup has no moving parts and is sensitive to weak sources 



(N 



(emitting single photons) of light as it relies on photon counting and has potential applications in both classical polarization measurements 
' and quantum communication scenarios. In our implementation, we are able to reconstruct the Stokes parameters of pure polarization 
^ states with an average fidelity of 99.9%. 
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Polarization measurement (also called polarimetry) is fundamental for many optical measurement tech- 

IT) ■ 

£^ ■ niques. In classical optics, applications include stress measurement, magnetic field sensing and optical thin 
q h ' film characterization [Q. More recently, polarization measurement of single photons prepared in particu- 

Ch ' lar ways has allowed experimental tests of non-locality and enabled applications in the field of quantum 

d \ 

q-i communication and information processing. 
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The purpose of this paper is to report on an implementation of an optimal polarimeter without moving 

X' 

' parts suitable for accurate polarization measurement at faint light (single photon) levels. Polarimetry 

_ _ i 

is usually performed using a combination of rotating wave plates and a linear polarizer. Implementations 
without moving parts to ensure speed and reliability have been reported before |2l3j and exhaustive research 
has been performed on optimizing them ^4 , 5JE1E1IH1 - Both types of polarimeters are applicable to complete 
state estimation techniques (also called state tomography) which are of interest in quantum information. 
Optimal polarimeters have been studied recently in the context of quantum key distribution |lUU12j where 
the concern is to obtain the best possible estimate on the polarization state of a limited collection of 
photons. While tomographic techniques are fairly common in quantum optics and related areas it 
should be noted that classical optics has used such techniques ever since Stokes showed how to determine 
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completely the polarization of an unknown state of light |13j . 

One instrumentally motivated way to fully characterize the polarization state of light uses three (Stokes) 
parameters S x , S y and S z . The three parameters may be augmented by an intensity S m and together form 
a Stokes vector, which can be normalized to be S = (1, S x /S m , S y /S m , S z /S m ); a reduced Stokes vector 
S r = (S x , S y , S z )/S m is used to identify a point in the Poincare sphere. 

A minimal scheme of measurement requires only four detector readings to estimate the Stokes param- 
eters. This corresponds to determining the overlap of the Stokes vector with four non-coplanar vectors 
in the Poincare Sphere. It has been shown that if this vector quartet defines a tetrahedron, then it has 
the necessary conditions for optimal polarimetry |IJ[71|S]. Independently, Rehacek et al showed that this is 
also true for two level quantum systems ^Uj. Optimal here is defined as having the least uncertainty for 
determining polarization from a given amount of light energy and that all possible input states are deter- 
mined equally well. In the context of photon counting, it means that the maximum amount of polarization 
information is extracted per photon counted. Thus, polarimeters based on a tetrahedron measurement are 
optimal and not only minimal. 

Independent of previous descriptions of four detector polarimeters, we designed and implemented a 
division-of-amplitude polarimeter based on the tetrahedron measurement. The setup is shown in Fig 1. This 
works by dividing the incoming light beam using a partially polarizing beam splitter (PPBS). Polarizing 
beam splitters are then used to analyze the daughter beams in different polarization bases. The design goal 
was to optimize the beam splitting ratio of the PPBS; optimizing this gives us an optimal polarimeter. 

Although previous reports describe division-of-amplitude polarimeters with fewer optical components |2j 
(no quartz compensators and no waveplates), they are often less than optimal tetrahedron measurements. 
Furthermore, these original designs place stringent demands on the performance of the PPBS. It must 
rotate light leaving the beamsplitter into the correct polarization basis as well as preserve the phase 
relationship between the horizontal (H) and vertical (V) polarizations. 

In our experiment we relaxed the requirement of having minimal optical components and aimed only to 
achieve optimal tetrahedron measurements without moving parts. We used a non-ideal PPBS, specifying 



Optimal Photon Polarimeter 
p PBS HWP at 22.5° 



Light of 

unknown 

polarization 

QWP at 45° 



Det 4,1 



/ Quartz 
Plates 



Det 1,1 
O 
Det 2, 1 



O Det3 ''3 

Figure 1. Setup for optimal polarimetry. The partially polarizing beam splitter (PPBS) light splitting ratio is important for our 
selection of waveplate angles. The quartz plates remove the relative phase difference between the horizontal and vertical polarizations 
introduced by the PPBS. In both transmitted and reflected arms, polarizing beam splitters then divide the light further between two 
detectors. Note that beams enter and leave the beamsplitter faces at 90° so alignment is straightforward. 

only an intensity splitting ratio required to work with selected waveplates (an ideal PPBS would not need 
waveplates and compensators but discussion of how to build such a device is out of the scope of this report; 
work continues in this direction j^j ) . Quartz compensators were used to correct for any relative phase shifts 
between the H and V polarizations introduced by the PPBS. Although a non-ideal PPBS requires us to 
have extra components they do not increase the complexity of the setup by much; in fact it is somewhat 
simpler because we are able to work with incident and exit angles of 90° for most of our elements. 

To understand how we found the optimal beamsplitting ratio, let us associate each detector with a vector 
bj from the tetrahedron quartet. Each vector bj is related to a measurement outcome Ij by 

j 

where Ij is the intensity of light falling on the associated detector j (for single photon counters, Ij refers 
to the photon count rate). The intensities may be expressed as a vector / = (Ii, I2, I3, 14), so we have 
I = B ■ S with an instrument matrix B that relates a Stokes vector, S to measured outcomes. When light 
of polarization bj is used, the normalized intensity of light falling on detector j is ^- At all other detectors, 
the intensity is g. Similarly, when light with polarization conjugate to bj is introduced, no light arrives at 
detector j but is split between all other detectors equally. 

The optimal splitting ratio for the PPBS is then most easily determined using Jones vectors (or the 
equivalent spin-i vectors for polarization qubits). Following Figure 2, suppose the PPBS has the property 
that input light with the polarization state (^) leads to the polarizations (^j) and (^) in the transmitted 
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and reflected arms, respectively. Two orthogonal vectors may be expressed as (J^^g) and ( _e cos s g me ) ■ 
Thus, the normalized intensity of light falling on detector j (e.g. j = 1,2) is 



hi it 



cos 9 e % ^ sin 9 



ay 
fix 



ay cos 9 + fixe sm.9 



h/It 



-aye 1 ^ sin 9 + fix cos 9 



(2) 
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Choosing a different measurement basis for detectors 3 and 4, we get: 



h/ It 



ax cos 9' + fiye sin 9' 



axe^' sin 9' + /% cos 9' 



(4) 
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Since the tetrahedron can be oriented arbitrarily, we choose for convenience 45° linear polarization 
(9 = 7r/4, 4> = 0) in the transmitted arm, and circular polarization (9 1 = 7r/4 0' = it/2) in the reflected 
arm. Using these values we find splitting ratios of 



2 1,1 2 1 1 

x = — 7=, y = F= 

2 2^/3 2 2^3 

The same result was shown independently by maximizing the instrument matrix determinant 0. Beam 
splitters with this nominal splitting ratio were obtained from a commercial optical coating facility. 

To use this polarimeter, the instrument matrix is determined by calibration. For this purpose, we pre- 
pared input light with four different polarizations whose Stokes vector are non-coplanar, giving a set of 
four simultaneous equations which are used to determine the elements in the matrix, B. This matrix takes 
into account losses and any non-ideal splitting ratio x, y of the PPBS. As a choice of calibration vectors, 
we select vectors that are maximally far apart in the Poincare sphere. The natural choice for this is again a 
vector quartet defining a tetrahedron. The normalized Stokes vector representation for our selected vector 
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quartet used to calibrate our polarimeters are 



5 



' 1 ^ 



/ \ 

1 



h 



V / 



zr. /2 



(7) 



In anticipation of future experiments, we prepared these states using light at 702 nm from a spontaneous 
parametric down conversion source [Tl]. The down converted light was collected into single mode optical 
fibres and has a spectral bandwidth of 4.74 ± 0.04 nm ^3]. We prepared the tetrahedron vector quartet 
using a Glan-Thomson polarizer (extinction ratio of 10 5 ) followed by a halfwave plate (HWP) and a 
quarterwave plate (QWP) allowing an arbitrary pure polarization state to be generated for calibration. 
Light detection is performed by passively quenched Silicon avalanche photodiodes. 

By this method the instrument matrix inverse was found to be 



/ \ 
1.015 0.988 0.981 1.013 



B- 1 



2.779 -3.156 0.231 0.302 



-0.138 0.149 -3.309 3.493 



2.004 1.458 -1.784 -1.865 



(8) 



The uncertainty in the elements has a maximum value of ± 0.007 deduced from propagated Poissonian 
counting statistics. Here we are limited only by the final accuracy of the state preparation apparatus. 

With the instrument matrix determined, it is possible to reconstruct the Stokes vector of any input 
state. This was tested by preparing a large number of pure states using our polarization state generator. 
For each input state, the Stokes vector was reconstructed from the detected outputs. We characterize the 
accuracy of the measurement by a fidelity F of the reconstructed state to the prepared state. The fidelity 
for two states is defined as F = {Tr[y'^7<^/5 rec ^/cr^]} 2 : , where a t h and p rec are the coherency matrices 
(or density matrices) of the theoretical and reconstructed states, respectively [TS|. For pure states, this 
reduces to \{S rec ■ Sth), i-e. the overlap of their normalized Stokes vectors. 
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Figure 2. Fidelity(.F) of reconstructed Stokes Vectors. The average fidelity is 99.9% while the minimum fidelity is 99.4%. Input states 
were generated from H polarized light transmitted through a half-wave plate and a quarter-wave plate. 



The map of this fidelity is shown in Figure 2. The average fidelity is 99.9% with a smallest value of 99.4 
± 0.4 % for circularly polarized light. This is attributed to imperfections in the waveplates of our state 
generator. It should be noted, however, that the range in fidelity can be accounted for by approximately 
4 mrad alignment error in the waveplates which is comparable with the specified accuracy of 5 mrad for 
the control motors. 

A careful inspection of Figure 2 shows that this optimal polarimeter has no large areas of low fidelity 
showing that it estimates all input states equally well. This is of particular importance in low light ap- 
plications with single photon sources. In our results, cumulative photon counts at each measured point 
was on the order of 10 5 (approximately 0.028 pJ). Further work to determine the asymptotic efficiency of 
tetrahedron polarimeters is currently planned. 

In conclusion, we have presented an experimental realisation of optimal polarimetry by photon counting. 
In the context of quantum information, this is a realisation of single qubit tomography based on the 
tetrahedron measurement (that is known to be optimal JO]). A pair of these polarimeters together could 
be used to perform state tomography of entangled photon pairs. 

The setup is very stable and fast as it has no moving parts. This could be applied to scenarios requiring 
time-resolved polarimetry or advanced quantum communication schemes, where optimal measurements on 
qubits can enhance communication bandwidth |17j . In its current form, the setup is able to reconstruct the 
Stokes vector of any input polarization to at least a fidelity of 99.4%. Although the work is presented for 
pure polarization states, there is no reason to believe it will not work for partially polarized light (called 
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mixed polarization states) to be addressed in future work. 
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